function main
im=double(imread('lena512.bmp'));
fim=fft2(im);
pcimg=img_polar(im,255,360);
fpcimg=img_polar(fim);

figure;
subplot(2,2,1); 
imagesc(im); colormap gray; axis image;
title('Input image');  

subplot(2,2,2);
imagesc(log(abs(fftshift(fim)+1)));  colormap gray; axis image;
title('FFT');

subplot(2,2,3);
imagesc(pcimg); colormap gray; axis image;
title('Polar Input image');  

subplot(2,2,4);
imagesc(log(abs(fpcimg)+1));  colormap gray; axis image;
title('Polar FFT');


end